Semiclassical spin liquid state of easy axis Kagome antiferromagnets 
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Motivated by recent experiments on Nd-langasite, we consider the effect of strong easy axis single- 
ion anisotropy D on S > 3/2 spins interacting with antiferromagnetic exchange J on the Kagome 
lattice. When T <g DS^, the coUinear low energy states selected by the anisotropy map on to config- 
urations of the classical Kagome lattice Ising antiferromagnet. However, the low temperature limit is 
quite different from the cooperative Ising paramagnet that obtains classically for T <^ JS^ . We find 
that sub-leading 0{J^S/D^) multi-spin interactions arising from the transverse quantum dynamics 
result in a crossover from an intermediate temperature classical cooperative Ising paramagnet to a 
semiclassical spin liquid with distinct short-ranged correlations for T ^ J S/D . 
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Introduction: The insulating magnet Nd-langasite has 
localized Nd'^"'" moments that live on sites of a 2- 
dimensional Kagome lattice made up of corner-sharing 
triangles (Fig [T^) . These moments carry a total angu- 
lar momentum Jion — 9/2 and are subject to a strong 
single-ion anisotropy term D ~ lOK [l|] that picks out 
the crystallographic c axis as the common easy axis of 
all the spins [3|- Although they interact with a size- 
able nearest neighbour antiferromagnetic exchange cou- 
pling J ~ 1.5K [l| (corresponding to a Curie- Weiss tem- 
perature Qcw = — 52K), Nd-langasite does not exhibit 
any magnetic order down to 50 mK [l|, IsJ, IJ] . Such be- 
haviour contrasts strikingly to conventional insulating 
solids with localized magnetic moments and concomitant 
short-ranged exchange interactions [5| — these usually en- 
ter a magnetically ordered state at low temperature, the 
nature of which can often be understood quite simply in 
terms of the classical energetics of the leading exchange 
interactions. 

When these interactions compete due to the geome- 
try of the lattice, as is the case in Nd-langasite, one of- 
ten obtains a large degeneracy of inequivalent classical 
ground states which prevents ordering and instead re- 
sults in unusual cooperative paramagnetic behaviour for 
a range of intermediate temperatures below Qcw - such 
systems are called frustrated magnets [6|. Frequently, 
quantum effects and sub-leading interactions eventually 
do lead to an (often complex) ordered state at still lower 
temperature — Kagome lattice magnets provide many ex- 
amples of this [Tl, Isl, Isl . In some other cases, such as the 
S = 1/2 Kagome magnet herbertsmithite [10|, there is 
apparently no tendency of the spins to form an ordered 
arrangement even at the lowest temperatures accessible 
to experiment. 

Systems such as Nd-langasite and Herbertsmithite pro- 
vide possible realizations of so-called spin-liquid states, 
which have been the subject of sustained theoretical ac- 
tivity pj| going back to the seminal work of Fazekas 
and Anderson [ij]. Much of this theoretical activity 
has focused on the challenging case of Heisenberg ex- 



change interactions (isotropic in spin space) and low spin 
S = 1/2 (in which quantum fluctuations are expected to 
be strongest) [11|. In the opposite Ising limit in which 
the exchange interactions only couple one component of 
neighbouring S = 1/2 moments on the Kagome lattice, 
there are no quantum effects and the system remains in 
a cooperative paramagnetic state with short ranged spin 
correlations all the way down to zero temperature [l3|. 

In order to model the case at hand, Nd-langasite, we 
consider moments with larger spin S > 3/2 interacting 
with isotropic Heisenberg exchange J on the Kagome lat- 
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tice, but which — unlike the case considered in 
subjected to a strong single-ion anisotropy D that picks 
out a common easy axis for all the moments: 
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When D dominates over J, the moments prefer collinear 
spin states that correspond to configurations of a clas- 
sical Ising antiferromagnet. However, quantum fluctua- 
tions induced by transverse components of the exchange 
coupling can lead to low temperature behaviour quite dif- 
ferent from the classical cooperative Ising paramagnet, 
and here we explore this possibility in some detail. 

Our results are readily stated: We find that such easy 
axis magnets do not develop any long-range order down 
to very low temperature, and are thus good examples 
of genuine spin-liquid behaviour in a system of quantum 
spins. More precisely, we find that there are two quali- 
tatively distinct regimes separated by a cross-over tem- 
perature T* - J^S/D'^. For T* <T <- JS'^, the short 
ranged spin correlations of the system are well-described 
by the correlations of the cooperative Ising paramagnet 
described above. Below T*, the leading effects of vir- 
tual quantum fiuctuations begin to dominate, leading to 
a qualitatively different semiclassical spin liquid regime 
in which the liquid structure factor of the spins encodes 
distinctive short-ranged correlations, but there is no long 
range order of any kind (this is in sharp contrast to 
easy axis Kagome antiferromagnets with 5* = 1 moments. 
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FIG. 1: (color online), a) A L^ ~ 4, Ly — 3 piece of the 
Kagome lattice, with principal directions To and Ti , and three 
sublattices of sites marked; also shown is a class of partially 
ordered states with macroscopic entropy that simultaneously 
minimize Ji and J2. b) Unfrustrated bonds of type n and 
mapping of Ising spins to dimers on the dual dice lattice. 



where the quantum dynamics connects different classical 
ground states and induces spin nematic order [l5|). 

Effective Hamiltonian: When D dominates over J, the 
system prefers collinear spin configurations with S^ = 
GiS (with (Ti — ±1). This degenerate ground-state man- 
ifold of the unperturbed problem can be thought of in 
terms of configurations of the Ising pseudo-spin variables 
a. The low energy physics in this regime is then best 
described in terms of an effective Hamiltonian Ti that 
encodes the splitting of this degenerate Ising subspace to 
each order in J / D. To order J^ /D^ , H is given as 16 | 
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where Ji = JS'^ , J2 = 4]j-f(2s-i)'^ ' ^^'^ ^^^ exchange field 
Hi = Tij aj with Fy = 1 for nearest neighbours and zero 
otherwise. In the above, the first term corresponds to the 
leading effect of the z component of the spin exchange, 
while the second term arises from virtual quantum tran- 
sitions of pairs of anti-aligned spins out of the low energy 
Ising subspace. An additional 0{J'^^ / D^^~^) pseudo- 
spin exchange term, representing real quantum transi- 
tions, is subleading for S > 3/2 [16|. 

Classical Ising regime: For temperatures T well above 
the exchange energy scale JS"^, but lower than the 
anisotropy energy DS^, the behaviour of the system will 
be that of to the high temperature paramagnetic regime 
of the classical Ising model on the Kagome lattice. As 
the temperature is lowered, the exchange begins to make 
itself felt, and the system crosses over to an intermedi- 
ate temperature regime J2 <C T ^ Ji whose physics is 
controlled by the ground states of the classical Kagome 
lattice Ising antiferromagnet (KIAF). 

In these ground states each triangle has exactly one 
frustrated bond (connecting a pair of aligned spins). This 
'minimally frustrated' ensemble of states has a residual 
entropy of 0.502^^ per site [13[, and the loss of en- 




FIG. 2; (color online). (a)Specific heat d, and uniform sus- 
ceptibility X in the crossover region, (b) The q^ = Qy cut for 
S{q) is also shown for two different system sizes {L — 48, 60) 
for two different temperatures (/3J2 ~ 0.0, 3.0). 



tropy coming from the high-temperature paramagnet is 
reflected in a featureless peak in the specific heat Cy (the 
area under the C^/T curve measuring this entropy loss). 

This ensemble of ground states can be conveniently 
represented by dimer coverings on the dual dice lattice 
wherein a dimer is placed on every link of the dual lattice 
that intersects a frustrated Kagome lattice bond (Fig[T]D); 
note that in this dimer model, a hard dimer constraint 
is operative only on the 3-coordinated dice lattice sites, 
while the 6-coordinated sites have a soft constraint that 
an even number of dimers touch them. 

As mentioned earlier, this ensemble [13| yields corre- 
lations that are extremely short-ranged and hence a fea- 
tureless spin structure factor (Fig[3K). The temperature 
regime J2 <C T -C Ji is thus a classical cooperative Ising 
paramagnet with a reduced (compared to the high tem- 
perature Ising paramagnet) value of the magnetization 
fluctuations as reflected in Tx (x being the linear sus- 
ceptibility to a magnetic field along the easy axis). 

Minim,izing J2: To understand the behaviour of the 
system at still lower temperatures, we first focus on the 
T = ground states of HeS and show that minimally 
frustrated configurations which also minimise the multi- 
spin interaction J2 satisfy the criterion that no spin be 
the minority spin of both triangles to which it belongs. 
To this end, we first note that the multispin interaction, 
projected to the minimally frustrated Ising subspace, as- 
signs a 'potential energy' —2J2{n— 1) to an unfrustrated 
bond of type n — here we classify an unfrustrated bond 
to be of type 3 (1) if the two spins it connects are both 
majority (minority) spins of the other triangles to which 
they belong, while an unfrustrated bond connecting one 
majority and one minority spin is of type 1 (Fig[T)3). 

As an unfrustrated type n bond of a minimally frus- 
trated Ising configuration is surrounded by a double 
rhombus of the dual dice lattice with n dimers on its 
perimeter (Fig[T)D), the projected multispin interaction 
J2 translates to a dice-lattice dimer model with interac- 
tions on double-rhombii surrounding unfrustrated bonds 
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where Nt is the number of triangles in the kagome lattice 
and /„ is the fraction of unfrustrated bonds of type n, 
i.e with n dimers on the perimeter of the corresponding 
double-plaquette (Fig[T}D.). 

Thus, up to a constant, the interaction energy counts 
the total number of dimers on the perimeters of all 
double-plaquettes corresponding to unfrustrated bonds. 
Since every single plaquette with n dimers on its perime- 
ter is part of 4 — n such double-plaquettes, the n dimers 
will contribute with multiplicity 4 — n to this total num- 
ber. Thus, if the fraction of single plaquettes with n 
dimers on their perimeter is 5„ (with J2n=o9n — 1j since 
each elementary plaquette can have 0, 1, or 2 dimers on 
its perimeter), the interaction energy of the configuration 
is —Np J2n=o "-(4~"-)3n, where Np = 3Nt/2 is the num- 
ber of elementary plaquettes. The interaction term thus 
assigns energies — quadratic in the number of dimers — to 
different types of elementary plaquettes of a dice lattice 
dimer configuration: 
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where \nP) denotes elementary plaquettes with n dimers 
on their perimeter. Since each dimer is on the perimeter 
of two elementary plaquettes, we have J2n=o ^9"- ~ ^/^j 
and this constraint along with X]n=o 9n = ^ allows us 
to minimize the potential energy Hp): Hp is minimized 
when go = (which fixes gi = 2/3 and 52 = 1/3) and 
this immediately gives the criterion that no spin be the 
minority spin of both the triangles to which it belongs. 

We have investigated the set of minimally frustrated 
configurations that satisfy this minority spin rule in some 
detail, and find that it is possible to construct a large 
subset of states (i.e. with macroscopic entropy) satis- 
fying this rule, and related to each other by local spin 
flips (Fig[T^). This construction immediately provides a 
lower bound of kp ln(2)/6 per site on the entropy of the 
ground states of Hp . The question then arises whether 
the correlations in the T ^ limit remain short-ranged, 
or whether order by disorder occurs here, and below we 
address this question numerically. 

Loop algorithm: In order to obtain reliable numerical 
results that can be used to settle this delicate question, 
we have developed a new loop algorithm that can handle 
the non-trivial Boltzmann weight associated with Hp as 
well as keep track of the different constraints on 3- and 
6-coordinated sites, while retaining the efficiency of the 
usual hard-core dimer model loop algorithms [iTJ, ll8| . 

The algorithm proceeds as follows: A change in dimer 
configuration is initiated by starting at a randomly cho- 
sen 6-coordinated site (say fio) and moving, with equal 
probability, to one of its six 3-coordinated neighbours, 
say i. If the link (fj-ai) is covered by a dimer, this dimer 
is rotated about the pivot site i, so that it now covers link 
(fii) — fi is chosen from the three possible 6-coordinated 
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FIG. 3: (color online). (a-b-c)Spin structure factor S{q) 
shown for three different temperatures, showing crossover to 
semiclassical spin liquid regime, d) Area under the CvjT 
curve gives the reduction in entropy associated with the 
crossover to the semiclassical spin liquid. 



neighbours of i according to a table of probabilities sat- 
isfying detailed balance (if ^ = ^0, the attempted up- 
date of the dimer configuration ends without making any 
change in the configuration). On the other hand, if the 
link (fJ-oi) is unoccupied by a dimer, and the dimer touch- 
ing i covers a different link (fii), this dimer is rotated 
about pivot i to cover link (fJ-oi) with a certain probabil- 
ity p chosen to satisfy detailed balance (conversely, with 
probability I ~ p, the attempted update ends without 
making any change in the dimer configuration). 

In both cases above, the first step of the update pro- 
cedure results in a violation of the (soft) constraint at /ig 
as well as at another 6-coordinated site /U. By repeat- 
ing this set of moves with probabilities chosen to satisfy 
detailed balance at each step, the site fj, can be moved 
around in a closed loop until it finally meets /iq again 
and heals the 'defects' that were originally introduced at 
that both sites. When this happens, a large change in 
the dimer configuration is effected along a closed loop of 
links, and this large change can be accepted with unit 
probability, thereby providing an efficient means of sam- 
pling the Gibbs distribution associated with Hp. 

The semiclassical spin liquid: From numerical simula- 
tions on Lx = Ly ~ L size systems using this algorithm 
(with L ranging from 10 to 60), we see no evidence at 
all of any phase transition as we lower the temperature 



to access the T — * limit. This is evident from the be- 
haviour of the specific heat per site, which converges very 
quickly with system size, and does not show any singular- 
ity in the thermodynamic limit (Fig[2^). In addition, the 
spin-spin correlators as well as the bond-energy correla- 
tors show no long range order at any wavevector down to 
the lowest temperatures we study (Fig [2b)- The system 
thus remains in a short-ranged ordered spin-liquid state 
down to the lowest temperatures. 

Although there is no phase transition, we find that the 
liquid state at low temperature is quite different from 
the intermediate temperature cooperative Ising param- 
agnet. This crossover to a distinct semiclassical spin liq- 
uid regime at low temperature (Fig [2K) is evident for 
instance in the temperature dependence of the specific 
heat per site Cy/Ng'- The Cu vs T curve shows a distinct 
but non-singular peak at T* k, I.3J2 that reflects the 
loss of entropy during this crossover from the cooperative 
Ising paramagnet to the low temperature limit in which 
the configurations sampled predominantly obey the mini- 
mum J2 constraint (from the area under the Cy/T curve 
(Fig [3Ji) and knowledge of the residual entropy of the 
cooperative Ising paramagnet, we estimate the residual 
entropy of the semiclassical spin liquid to be 0.32fcB). 

A clear signature of this crossover to the semiclassical 
spin liquid regime below T* can be obtained by monitor- 
ing the spin structure factor (that can be probed in neu- 
tron scattering experiments) S{q) — \So{q)ex.p{iqy/2) + 
Si{q) + S2{q)exp{iqx/2)\'^ where Sa{q) is the Fourier 
transform of the spin density on sublattice a of the 
Kagome lattice and qx (%) refers to the projection of 
q on to lattice direction Tq (Ti) measured in units of in- 
verse Bravais lattice spacing (Fig[T^). From Fig [3] a), b), 
c), we see that the structure factor evolves continuously 
from being quite featureless in the classical cooperative 
Ising regime T* <C T <C Ji to developing characteristic 
crescents of high intensity diffuse scattering in the low 
temperature semiclassical spin liquid regime T <^ T*, 
with precursors of these features being already present 
at T ~ 2T*. In addition, this crossover is also charac- 
terized by a change in the magnetization ffuctuations as 
reflected in the value of Tx (Fig [5^). 

Experiments on Nd-langasite: As mentioned earlier, 
recent experiments on the spin-9/2 easy axis Kagome an- 
tiferromagnet Nd-Langasite have seen a liquid-like state 
with fluctuating moments and no long range order down 
to bOmK. From the estimated values ^1]] of the isotropic 
exchange interaction and the single-ion anisotropy (J ^ 
1.5K, D ~ WK), the crossover from a classical coop- 
erative Ising paramagnet to a semiclassical spin liquid 
only occurs for temperatures significantly below 50mK 
(T* « 16mK) in this model with an isotropic exchange 
and a simple single-ion anisotropy term (which is ex- 
pected [l| to be a good starting point for Nd-Langasite) . 

We therefore expect a simple classical Ising description 
to work fairly well for the bulk of the temperature range 



studied in the recent low temperature experiments. Since 
correlations in the cooperative Ising paramagnet are ex- 
tremely short-ranged and featureless, this is consistent 
with the fact that data from recent neutron scattering 
experiments (which probe the spin structure factor) show 
nearly featureless diffuse scattering that can be fit quite 
well to a model [SJ, |j] of spin correlations in which nearest 
neighbour spins are correlated, but there are no correla- 
tions of spins further away from each other. 

We hope that our results provide further motivation to 
study this magnet at still lower temperatures, at which 
the crossover to the interesting low temperature semi- 
classical spin liquid should become apparent. Another 
possible avenue for exploring this crossover in greater de- 
tail involves identification of other easy axis Kagome an- 
tiferromagnets in which the separation of scales between 
J and D is not so large (or the overall scale of both J 
and D is somewhat larger) so that the crossover to the 
semiclassical spin liquid occurs at more easily accessible 
temperatures and we hope that the results of our work 
provide motivation for exploring this possibility as well. 
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